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Abstract 

The consistency of quantum adiabatic theorem has been doubted re- 
cently. It is shown in the present paper, that the difference between the 
adiabatic solution and the exact solution to the Schrodinger equation with 
a slowly changing driving Hamiltonian is small; while the difference be- 
tween their time derivatives is not small. This explains why substituting 
the adiabatic solution back into Schrodinger equation leads to 'inconsis- 
tency' of the adiabatic theorem. Physics is determined completely by 
the state vector, and not by its time derivative. Therefore the quantum 
adiabatic theorem is physically correct. 

0.1 Introduction 

Quantum adiabatic theorem (QAT) dates back to the early years of quan- 
tum mechanics[l]. It has important applications within and beyond quantum 
physics. In 1984, M. Berry found there is a geometrical phase in the adiabati- 
cally evolving wavefunction besides the dynamic phase[2]. B. Simon pointed out 
Berry's phase factor is the holonomy of a Hermitian line bundle [3]. This started 
a rush for geometrical phases in quantum physics [4], which helped people to 
get deeper insight into many physical phenomena, such as Bohm-Aharanov ef- 
fect, quantum Hall effect, etc. Recently, the quantum adiabatic theorem has 
renewed its importance in the context of quantum control and quantum compu- 
tation^]. More recently, however, the consistency of QAT has been doubted[6]. 
In their paper entitled "Inconsistency in the application of the adiabatic the- 
orem", K.-P. Marzlin and B.C. Sanders gave a proof of inconsistency of QAT, 
and declared that the standard treatment of QAT alone does not ensure that a 
formal application of it result in correct results. This interesting suggestion has 
attracted attention from the physics circle [7]. The purpose of this letter is to 
point out that the QAT does give approximate state- vectors (wavefunctions), 
but not necessarily the approximate time derivatives of state-vectors. While 
physics is completely determind by the state- vector, and has nothing to do with 
its time derivative. Therefore the QAT is physically correct. What leads to 
'inconsistency' of the QAT is neglecting the fact that the adiabatic approxi- 
mate state-vector does not necessarily give the approximate time derivative of 
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state- vector (\\ip A (t) - ip(t)\\ < 1 ^> 
denotes the norm of state- vector ip ). 



~ V'W < 1, where ||<p|| = VT^) 
0.2 Standard treatment of QAT 

Suppose that the Hamiltonian depends on N real parameters R 1 , . . . , R N : 

H = H(R 1 , . . . , R N ) = H(R) (1) 

When the representing point of the Hamiltonian describes slowly a finite curve 
C on the N— dimenssional parameter manifold M 

C : R a = R a (t), VtE [0,T],1 < <7 < TV (2) 

where T is the evolution time, let us study the evolution of the system. The 
instantaneous Hamiltonian's eigen equation is 

H(R)u n {R) = E n (R)u n {R) (3) 

Getting to the rotating representation 

m = E Cn(t)u n (R(t)) cxp[-^ f E n {R(t'))dt'] (4) 

n>0 h J ° 

we get the Schrodingcr equation 

Cm{t) = -E„>0 (u m (R(t)) | U n (R(t))) X 



cxp |l jf - ^„(i?(t'))]^'| c„(i) (5) 

To avoid confusion of infinitesimals of different orders and to show what 'rapidly 
oscillating' means, let's change to the dimensionless time r = t/T 
£c m (r) = - £„>o (u m (R(T)) | £u n (R(r))) x 



exp jir £[E m (R(r')) - B„(A(r , ))]dr , | c„(r) 



(6) 



where 

c m (r) = c m (Tr) = c m (i), %) - fl(Tr) = i?(t) (7) 

The initial value problem of the above differential equations is equivalent to the 
following integral equations 

c m (r) = c m (0) - £„> /; I £u n (R( Tl ))) x 



exp|lT^ Tl [£^(A(T , ))-£; n (A(T , ))]dT , |g n (Ti)dTi 



(8) 
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Let's slow down evenly the changing speed of the Hamiltonian while keep the 
finite curve C fixed. Mathematically, that is to let T — * oo, while keeping the 
function form of R{t) unchanged. The oscillating factors in the integrand ensure 
vanishing of the corresponding integrals. There is no resonance problem in the 
mathematical context. For the practical physical problem, slowly changing of 
the Hamiltonian means T is such a long time that 



(R(t)) \u n (R(t)))h 



E m (R(t)) - E n (R(t)) 



« f (9) 



In both mathematics and physics contexts, the integral equation (8) can be 
approximately rewritten as 

ci(r) = c m (0) - J (u m (R( Tl )) | ^-u m (R( Tl ))j ^(Ti)dn (10) 
Solving this equation by using iteration gives 

c^(r) = exp|-^ (u m (R( Tl )) \ ^u^RiT,))^ dr^ c m (0) (11) 
This proves the QAT. 

0.3 Analysis of 'Inconsistency' of QAT 

When we substitute the adiabatic approximate solution (11) back into the in- 
tegral equations (8), the equations approximately hold. 

« - Hn^ra) Jo I £«n(%l))> >< 

exp|^T^ Tl [i; ro (%'))-^n(%'))]^'}^(Ti)rfri (12) 

However, when we substitute the adiabatic approximate solution (11) back into 
the differential equations (6) whose initial value problem is equivalent to the 
integral equations (8), we obtain 

»-En(H <«m(%)) I X 

C ^\^T J\ Em {R(T')) -E n (R(r'))]dr'^c^r) (13) 
Considering that V(0) can be an arbitrary state- vector, we have 

w (u m (R(T)) | ^-u n {R{r))^ ,Vm^n (14) 

which is false. Notice that the right-hand side of (13) is the derivative of the 
right-hand side of (12), while (12) is correct and (13) is incorrect. In order to 
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understand the situation we are facing, let's study the following basic mathe- 
matical fact. Let \ip(t)) = \0)er lut +e\l)e- lt / {e2 \($ < e < 1), \<p(t)) = \0)e~ ilvt , 
where |0), |1) are eigen vectors of the 1-dimensional harmonic oscillator energy. 



||^))-|^))||=s«l, .:W))K\<P(t)) 



While 



= i/ £ »i,.-.iv(i))^i^)) 



(15) 
(16) 



The above example shows that two approximately equal time-dependent state- 
vectors do not necessarily have approximately equal time derivatives. Therefore 
the approximate solution to integral equations (8) does not ensure that the 
equivalent differential equations (6) approximately hold. It is neglect of this 
basic mathematical fact that leads to 'Inconsistency' of QAT in [6]. 

QAT gives the approximate state-vector, not the approximate time deriva- 
tive of the state-vector. All the physics is, however, determined by the state- 
vector itself, not by its time derivative. Therefore QAT is completely correct 
physically. 

0.4 An exactly solvable example 

Let's consider an exactly solvable example, the evolution of the spin wavefunc- 
tion of an electron in a slowly rotating magnetic field B{t) — B a ( i cos + 
j sin %p-t). The instantaneous Hamiltonian is 

H(t) = -7? • B(t) = ±? ■ B(t) = ||V • B(t) 



ehB 
2m 



e 

i2nt/T 



i2nt/T 





e 

i2irt/T 



i2-Kt/T 





Its eigenvalues are E± (t) — ±e. And the corresponding eigenvectors are 



u±{t) = T2 



-iirt/T 



± e mt/T 

The exact general solution to the Schrodingcr equation 

d 



(17) 



(18) 



(19) 



or 



is 



ih— 
dt 









= e 



x(t) 

y(t) 



e 

e i2nt/T o 



27rt/T 



x(t) 

y(t) 



-A~ 1 [a(B + CY ct +c 2 (B- C)e~ lCt ]e 



[ Cl e iCt + c 2 e- iCt ]e 



(20) 



(21) 
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where A = e/h,B = tt/T, C = V A 2 + B 2 , and c\, C2 are the integral constants. 
The specific solution determined by the initial condition 



V(o) 



" x(0) ' 


1 


' i " 


. y(o) . 


= 71 


i 



IS 



V2 



(cos Ct-i^sin Ct) 
(cos Ct-i^sinCi) e" 



Let's get into the rotating representation. 

iP(t) = c+(t)u + (t)e- tAt + C-(t)u_(t)e lAt 

The exact Schrodinger equation becomes 

c+(t) = 8 Be !2if c_(t) 



c_(*) =*£e M*) 



Its general solution is 



" c+(t) ' 









£ iAt ( C-A c r c iCt 



- c + j A c"c~ iCt ) 



The specific solution determined by the initial condition 



" MO) " 




' i " 


. MO) . 








IS 



(cos Ct 



i£ sin Ct) e iAt 



(igsin Ct)e 



iAt 



(22) 

(23) 
(24) 

(25) 
(26) 

(27) 
(28) 



The adiabatic approximation means neglecting the non-diagonal (n ^ m) terms, 
which contain oscillating factors, on the right-hand side of differential equations 
(25). The adiabatic approximate solution detcrmind by the initial condition 
(27) is 

! (29) 



" 4(t) - 




' 1 ' 


c A (t) 








Getting to the dimenssionless time r = t/T, we rewrite (28) and (29) as 
Mr) 



c_(r) 
y/(eT/h) 2 + tt 2 t - 



cos 



eT/h 



sin J (sT/h) 2 +it 2 t e ieTT l n 



S— sin y/ieT/hf + n^ e^/H 



(30) 
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It's easy to see that 





cf(r) 




" 1 






gf(r) 









" Mr) 




1 " 




cf(r) " 


. Mr) 












(31) 



(32) 



The difference between (30) and (31) is small, but rapidly oscillates with the 
dimcnsionless time r. Therefore, it's to be expected that the derivative with 

r of the difference is no longer small. In fact, letting F = \J (sT/H) 2 + it 2 , we 
have 



dT C - 



-to 



ireT/h 
F 



sin Ft) e isTT l h 
sin Ft + in cos Ft) e - ieTT / h 





i7re -i2eTr/h 



ci(r) 
'to 



(33) 



0.5 conclusion 



The above discussion shows: (i) The adiabatic state- vector ip A (t) does not satis- 
fies approximately the Schrodinger differential equation, but it satisfies approx- 
imately the equivalent integral equation, (ii) The QAT is completely correct 



physically. This is ensured by 



<C 1. But it's not necessarily 



true that 



< 1. (iii) Taking ip (t) for ip(t) will probably lead to 



contradiction. 

Even though we don't agree with [6], we still think it's an interesting work. 
Because it has raised an important question: In theoretical reasoning, one has 
to bear in mind that approximately equal functions do not have to have approx- 
imately equal derivatives. 

Acknowledgement 1 We thank L.-A. Wu for letting us know the work [6] 
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